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Abstract: We discuss non-perturbative QCD contributions to jet observables, computing 
their dependence on the jet radius R, and on the colour and transverse momentum of the 
parton initiating the jet. We show, using analytic QCD models of power corrections as 
well as Monte Carlo simulations, that hadronisation corrections grow at small values of 
R, behaving as 1/R, while underlying event contributions grow with the jet area as R 2 . 
We highlight the connection between hadronisation corrections to jets and those for event 
shapes in e + e~ and DIS; we note the limited dependence of our results on the choice of 
jet algorithm; finally, we propose several measurements in the context of which to test or 
implement our predictions. The results presented here reinforce the motivation for the use 
of a range of R values, as well as a plurality of infrared-safe jet algorithms, in precision jet 
studies at hadron colliders. 
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1. Introduction 

Jets play a key role in many experimental studies at the Tevatron, and will continue to do 
so at the LHC. They are useful, for example, in top reconstruction and mass measurements, 
in searches for the Higgs boson and new physics signals, and of course they are instrumental 
for QCD studies, such as inclusive-jet measurements, which in turn are an important input 
for the determination of parton distribution functions. 

Jets, however, are fundamentally ambiguous objects, reflecting the fact that the di- 
vergences of QCD perturbation theory make it impossible to define clusters of final state 
hadrons uniquely assigned to individual hard partons emitted at leading order. Most pro- 
cedures to cluster measured hadrons into jets, i.e. jet algorithms, deal with this ambiguity 
by introducing a resolution parameter to define when exactly a jet with substructure should 
actually be resolved into two jets. For currently used hadron-collider jet algorithms the 
main parameter is generally a threshold on the allowed opening angle of the jet, and is 
called R or jet radius (defined on the azimuth-rapidity cylinder). 

Given, on the experimental side, the increased precision required for present hadron 
collider studies, and the complex hadronic environment expected in LHC final states, and 
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given many years of theoretical advances in our understanding of QCD, it is both useful 
and possible to work towards a better quantitative understanding of how the choice of 
the jet algorithm and the specification of its parameters influence the clustering of QCD 
radiation, and our general understanding of the final state in hadron collisions. 

One can consider three classes of QCD effects: perturbative radiation, hadronisation 
and the underlying event. It should be realised that there is no unique, fully consistent, 
gauge-invariant way to distinguish these three classes: rather, one must try to define per- 
turbation theory, and the underlying factorisation scheme, precisely enough, in order to 
identify non-perturbative hadronisation corrections without double counting; underlying 
event corrections will, at some level, inevitably mix with hadronisation: one must argue 
that this mixing does not occur, or is parametrically controllable, at least at leading and 
next-to-leading power in the hard scale. In this respect, as we will see, the radius depen- 
dence is a very useful tool. With this premise, one would like to understand how each of 
these classes of contributions affects the momentum of a jet, as a function of the jet radius 
R, as a function of the centre-of-mass energy, and as a function of the properties of the 
parton that initiated the jet, such as its colour charge and its transverse momentum. 

These issues can be studied through parton-shower Monte Carlo generators, taking 
advantage of their detailed and well tested non-perturbative modelling. It can, however, 
be difficult to extract simple analytical understanding from such models. An alternative 
approach is to carry out analytical calculations directly. Perturbatively, one may work 
with soft and collinear approximations for gluon emission, and perform a resummation. 
The tools for this task, in the intricate case of jet production in hadron-hadron collisions, 
have been available for a long time [1-3], and phenomenological studies of the effects of soft 
and collinear logarithms near threshold have been performed, both for inclusive jet cross 
sections [4,5] and for jet shapes [6]. It is well understood that these resummations carry 
nontrivial information on the parametric size of non-perturbative corrections, as discussed 
in [7]. Furthermore, one can carry out non-perturbative studies using renormalon and 
related approaches [8-11]; these suggest that the dynamics of Monte Carlo models can in 
some cases be reduced, for sufficiently inclusive observables, to a single non-perturbative 
parameter, multiplied by analytically calculable coefficients. These models have had con- 
siderable success in the context of e + e" and DIS (for reviews, see [12,13]); in contrast, for 
QCD final-state observables at hadron colliders, only a handful of results exist, notably for 
the jet energy-flow profile [6], the out-of-plane momentum for Z+jet production [14] and 
away-from-jet energy flow [15]. 

It is perhaps worth emphasising that the non-perturbative corrections we are dis- 
cussing, especially at moderate^, have an impact that is comparable to that of higher-order 
perturbative effects, and may affect in a crucial way precision studies for many observables, 
even in the nominally asymptotic energy regime probed by the LHC. This is a familiar issue 
from QCD studies at LEP and HERA, especially those concerning event shapes. In the 
context of hadron collider phenomenology, non-perturbative radiation from hadronisation 
and underlying event is, for example, the main source of uncertainty in the determina- 
tion of the jet energy scale. This plays a key role in precision studies, as was shown at 
Tevatron [16]: each percentage point of uncertainty in the jet energy scale translates, for 
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example, into a 1 GeV uncertainty in the top mass, and into a 10% uncertainty in the 
single-inclusive jet pt distribution at p t ~ 500 GeV. In general, the steep slope of jet distri- 
butions amplifies the effects of comparatively small energy shifts due to non-perturbative 
radiation, and propagates their impact all the way to the highest energies [17]: this effect 
will actually be magnified at LHC, where the underlying event is expected to provide an 
energy density much higher than that measured at Tevatron. 

In this paper, we begin to address quantitatively a number of issues related to non- 
perturbative corrections to jet observables at hadron colliders. We focus mostly on the 
radius dependence of non-perturbative effects, but examine also the dependence on the 
quantum numbers of the partons initiating the jets, and the effects of the choice of jet 
algorithm, thus complementing and extending the more qualitative statements found for 
example in [18, 19]. For the sake of completeness, and in order to be able to compare 
perturbative and non-perturbative effects as precisely as possible, we begin in Section 2 by 
reviewing and organising perturbative results at NLO. Similar results have been derived 
and used before in order, for example, to allow for a precise matching between resummed 
and NLO jet cross sections [5]. In Section 3, we move to the main subject of our analysis, 
and provide an analytical estimate for the change Spt in the transverse momentum of a jet 
due to hadronisation, at leading power in pt- We employ the techniques of the dispersive 
approach to power corrections, adapted to the environment of hadron collisions. Our main 
result concerns the radius dependence of Spt, which behaves like 1/R at small R, in sharp 
contrast to the behaviour of underlying event contributions, which are expected to grow 
with the jet area and become negligible for very narrow jets. The 1/R dependence of Spt 
could have been deduced from the analysis of [7], predicting a 1/(QR) correction to the jet 
cross section. This fact seems however not to have been widely appreciated. Relative to [7], 
here we work with fully-defined jet algorithms, we show that corrections of O (1) vanish, 
while those of O (R) have a small coefficient, we provide the relation between the coefficient 
of 1/R and the non-perturbative parameter ao used in event shape studies (cf. [13] and 
references therein), separately for quark- and gluon-induced jets, and we outline the class 
of jet algorithms for which this relation is expected to be exact. 

In Section 4 we compare our analytical results with Monte Carlo simulations, while 
varying R, the main parameter entering the jet definition. We employ Pythia [20], as well 
as Herwig [21] (with Jimmy [22] providing the underlying event) to generate hadron level 
jet events, and we reconstruct the jets using different (IR safe) algorithms, for different 
parton channels and different centre of mass energies. It is possible to generate separately 
the hadronisation and underlying event components of Spt with Monte Carlo simulations; 
the main result for the hadronisation contribution is a striking confirmation of our analytic 
expression: the four different jet algorithms that we employ all agree with each other and 
with the analytic result, both in shape and normalisation, within margins that are quite 
reasonable given the inherent approximations of both models. The underlying event com- 
ponent is certainly much less well understood; it is clear, in any case, that the momentum 
shift due to the underlying event grows with R, and for moderate values of R it is several 
times larger at LHC than at the Tevatron, as might be expected. 

We conclude, in Section 5, by outlining some experimental consequences of our results 
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and discussing some possible measurements that could validate them. Specifically, we 
discuss the optimal choice of the jet radius R in different experimental circumstances, a 
practical way to measure the pt shift due to hadronisation and underlying event, and the 
effect of such a shift on the single-inclusive jet pt distribution. 

We regard these results as a first step towards a detailed analysis of jet physics beyond 
fixed order perturbation theory, which could be improved upon in several ways: our analytic 
estimate could for example be made more precise by matching it to a resummed prediction 
for a specific jet distribution; in general, we expect that beyond the leading 1/R behaviour 
non-perturbative radiation will not simply shift the cross section but also change its shape, 
as seen for example in [15]. One also expects that the shift will not be universal across 
different observables and algorithms, since typical jet algorithms introduce nonlinearities 
that are likely to spoil the simple pattern of exponentiation that underlies the power of 
one-gluon results [23]: the effects of this breaking of universality could be studied both 
with analytical and Monte Carlo tools. 

We believe, in any case, that our results strongly suggest that experimental collabo- 
rations should try to maximise the flexibility of their choices concerning jet analyses. A 
wealth of information, and perhaps even opportunities for discovery, can be missed, if fu- 
ture jet studies are confined to just one or two algorithms, and a handful of parameter sets, 
chosen ahead of time. LHC is a discovery machine, and in order to fully exploit its potential 
we must be prepared with a range of flexible tools, capable of meeting unexpected, as well 
as expected challenges. 

2. Perturbative R dependence 

We begin by reviewing the perturbative dependence on the jet radius R, for small val- 
ues of R, where a priori one can expect a logarithmic enhancement originating from the 
collinear singularity that is approached when the jet becomes very narrow. Specifically, we 
observe that perturbative lni? terms (discussed on many occasions previously, [5,24-27]) 
originate from the same phenomenon as the 1/R growth that we will find for hadronisation 
contributions: as the jet becomes narrow, partons radiated outside of it (in other words 
partons which are not recombined with the jet by the chosen jet algorithm) are allowed 
to become more and more collinear to the emitter, approaching the collinear-singular con- 
figuration. Since these contributions are uniquely associated with the outgoing jet, and 
are independent of the other hard emitters, they can be simply computed in the collinear 
approximation. 

In order to perform explicit calculations, we need to pick specific observables. Let 
us first consider the loss of transverse momentum for a leading jet, which is of relevance 
for example when using the hardest pair of jets in an event to reconstruct the mass of a 
decayed massive particle. As an example, consider the quasi-collinear branching of a quark 
with transverse momentum p t , which splits into a quark carrying a fraction z of the initial 
momentum, plus a gluon carrying the remaining fraction 1 — z. Given that the probability 
of such a quasi-collinear branching is just the corresponding splitting function, one can 
easily write down an expression for the leading perturbative contribution, at small R, to 
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the average change in transverse momentum of a high p t jet. One finds 
(Spt) P crt=Pt J^-J dz{max[z,l-z]-l) <° P qq ( z )Q(6- f aig (z)R) . (2.1) 

Here we constructed the change in transverse momentum as the difference between the pt 
of the leading jet (which can be a quark or a gluon jet, depending on whether z or 1 — z is 
larger) and that of the initial quark, accounting for the splitting probability P qq {z), averaged 
over the essentially collinear branching. The 6 function constraint denotes the condition 
for non-recombination of the softer parton into the leading jet, and / a i g (-z) is a function that 
depends on the jet algorithm. In the kt [28,29] and Cambridge/ Aachen [30] algorithms (for 
detailed definitions, see section 4), for example, one merely requires that the small angle 9 
between the quark and gluon be greater than R, so that fk t ( z ) = fcam(z) = 1. For stable- 
cone based algorithms, such as SISCone [31] 1 , one is required, in order to enforce stability, to 
first construct the energy-weighted centroid of the quark-gluon system, nj = zn q + {l — z)n g , 
where n q ^ are unit vectors along the the quark and gluon directions. One must then ensure 
that the angle between tij and the softer parton be greater than R, which implies that the 
parton will not be recombined. Projecting, for example, the unit vector along the jet 
hj = nj/\ftj\ onto the gluon direction, in the small-angle limit one finds 

6 jg = z8 qg , (2.2) 

relating the opening angle between the gluon and jet directions, 0j g , to the quark-gluon 
opening angle 6 qg = 6. Non-recombination, for the case where the gluon is the softer 
parton, so that z > 1/2, requires 6j g > R, which in turn gives 

9 qg >- z = R (l + ±^) ; (2.3) 

if on the other hand the quark is softer (z < 1/2), the corresponding result is obtained by 
replacing z with 1 — z in Eq. (2.3). This leads to 

fconc(z) = 1 + min (j-r^, ^~r) ■ (2 ' 4) 

As a consequence, for cone algorithms, the 6 integral in Eq. (2.1) is cut off at R/z (R/ (1— z)) 
for z > 1/2 [z < 1/2). The z-dependence of the integration region does not affect the 
coefficient of mi?, but gives an ^-independent shift, computed below. Note that this 
result is independent of the overlap threshold in the split-merge stage of SISCone and 
similar algorithms. 

Let us now return to Eq. (2.1). The logarithmic behaviour comes from the 9 integral, 
which is cut off by the O function at small angle, while the upper limit is given by the 
large-angle hard partonic structure of the event, and does not affect the R dependence in 
the present approximation. For a quark-initiated jet, define then 



L q = J dz min (z, 1 - z) P qq {z) = C F ^21n2 - ^ , 



(2.5) 



1 Kt this order the Midpoint algorithm [32] behaves in the same way, however it sutlers from infrared 
unsafety at higher orders. 
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while for a gluon-initiated jet the corresponding quantity is 

Lg = J dzmm(z,l-z) (±P gg (z) + n f P qg ^) = C A (2 In 2 - + n f T R L , (2.6) 

where we have included a 1/2 factor for Bose symmetry in the gg channel, and Pij{z) are 
the real emission parts of the leading order DGLAP splitting kernels. It is easy to verify 
that, for a jet originated by a parton of type i, Eq. (2.1) yields 



(SptYpert _ «£ 

Pt K 



Li\nR + 0(a s ), (2.7) 



in a fixed-coupling approximation, with corrections which are non-singular as R — > 0. The 
main feature of Eq. (2.7) is the logarithmic dependence on the radius R, which may spoil 
the convergence of the perturbative result at small R. For a complete description of the 
perturbative result at the smallest R values, one might need to resum this logarithmic 
enhancement to all orders, an analysis which we postpone to future work. 

For completeness, we also compute the ^-independent shift in transverse momentum 
between the stable-cone-type and kt (or Cambridge) algorithms. It is defined by 

(^r C )pcrt-(^)pcrt a, (2g) 
Pt VT 

where again i labels the parton species and we find 



K q = J ' dz min (z, 1 - z) In (f CO ne(z)) P qq (z) 



— + - In 2 + In 2 2 J C F ~ 0.323C F , (2.9) 
16 8 J 

for quarks, while for gluons one has 

K g = J dzmm(z, 1 - z) In (/ co „ c (^)) + n f P qg ) (2.10) 

1321 133 \ /241 25 \ 

1152 + "96" ln 2 + ln ' 2 J ° A + { 576 " 48 ln 2 J ^ ~ + 0m7n f TR " 

Numerically, one notes that the Ki ~ 0.3Lj. This is the cause of the feature (origi- 
nally observed in [29]) that, perturbatively, kt and cone algorithms behave similarly when 
\nR kt ~ 0.3 + lni? 

cono ot equivalently Rk t — l-35i? cone . 
A similar analysis can be carried out for other observables as well. One can for example 
consider the single-inclusive jet p t distribution. Assuming a spectrum that falls off as 1/p™ 
at Born level, one finds at NLO a ln R correction (in this case multiplying the Born level 
distribution oc l/pf) which is still of the form of (2.7), but with new coefficients L q (n) and 
L g (n). Defining as usual 

^) = ^mr(n) = - 7B + X;l (2.11) 

p=i 
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one finds 



L, 



, q (n) = [ dz (z^ 1 + (1 - z) n - 1 - 1) P qq {z) 



Jo 



= -C F [ 2iP(n - 1) + 2 7jB + 




) 



(2.12) 



as well as 




11 



6 



) 



(2.13) 



This pattern is general: for a generic jet observable, the logarithmic behaviour in R is 
dictated by collinear dynamics, and can be computed at NLO (and beyond) using only 
splitting function information rather than the full squared matrix element. In this regard, 
we note that Ref. [19] has discussed the R dependence of the average squared jet mass, 
observing that to first order in a s it behaves as a s R 2 p\. 

3. An analytical estimate for hadronisation corrections 

To be concrete, let us choose a definite process and observable. We consider single-inclusive 
jet production near partonic threshold in hadronic collisions, a process of current interest at 
the Tevatron and soon at the LHC. Our definition of threshold is that the scaled transverse 
momentum xt = 2pt/^/s approaches its kinematical limit x% = 1, so that all the available 
collision energy is converted into the jet transverse momentum. We have in mind the 
rapidity-integrated distribution, which however is known to be well approximated near 
threshold by the value of the distribution at vanishing rapidity [5]. We work then in the 
partonic centre-of-mass frame, and place the trigger jet at vanishing rapidity, with no loss 
of generality. At Born level we have the kinematics 




(3-1) 



where pj is the four-momentum of the trigger jet, pt its transverse momentum and p r the 
four-momentum of the recoil jet. At Born level both jets correspond to massless particles 
and it follows that 



Consider now the change in pt induced by a soft gluon emission. We must separately 
consider two alternative scenarios: either the soft emission is recombined by the jet al- 
gorithm with the hard parton associated with the trigger jet, or it is left unrecombined, 
so that the measured jet remains massless at this order. The recombination scheme we 
choose is the one commonly used for jet studies at the Tevatron, the E scheme, where 
the four-momentum of the jet is obtained by adding the four-momenta of the constituent 
partons. 

In the case where recombination happens the final state kinematics becomes 

Pi = (^P 2 t+Mf,p t ,0,0) (3.3) 

Pt = {Pt,~Pt, 0,0) , 

where the trigger jet now has a mass Mj, while the recoil jet is still a massless hard parton. 
Using energy conservation one has then 



pt + ^p 2 t + M] = ^s. (3.4) 

Expanding to first order in the mass Mj (since the gluon we recombine with the jet has 
been assumed to be soft) we obtain 

The change in p t from its Born value is given by 

r +/ M M j P3-k 

^ + (*) = -Vi = — (3 - 6) 

where the + superfix denotes the case when the gluon is recombined with the jet, and we 
wrote the jet mass as Mj = (p% + k) 2 = 2^3 • k, with p% the four-momentum of the hard 
massless parton initiating the jet, and k that of the soft gluon. 

When the gluon is not recombined, the recoil system is massive and one similarly finds 

S Pt (k) = ~ Ws = , (3-7) 

where P4 denotes the four-momentum of the massless hard parton recoiling against the 
trigger jet. Parametrising the soft gluon four-momentum as 

k^ = k t (cosh rj, cos 4>, sin <f>, sinh 77) , (3-8) 

where kt, rj and <f> are respectively the transverse momentum, rapidity and azimuth defined 
with respect to the beam direction, it is easy to see that to first order in the small quantity 
kt one has 

kt 

5p t (k) = — — (cosh 77 =F cos 4>) . (3-9) 

To obtain the average change in jet pt due to soft emissions, one has to multiply the result 
in Eq. (3.9), for the change in p t due to a soft emission, times the probability of emitting 
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the soft gluon, and subsequently integrate over phase space. In the eikonal approximation, 
appropriate to soft emissions, the probability (squared matrix element) for the emission 
of a gluon with momentum k from an ensemble of hard partons with momenta pi may be 
expressed as a sum over contributions from all possible colour dipoles 

\M\ 2 = \M \ 2 ^2d J W lj (k) , (3.10) 

where the sum runs over all distinct pairs (ij) of hard partons, or equivalently, as stated 
before, over all dipoles. The quantity |A4o| 2 is the squared matrix element for the hard 
scattering, which in our case has to be computed for each separate partonic subprocess 
contributing to the jet distribution, and contains the dependence on parton distribution 
functions. The contribution of each dipole Wij is weighted by the colour factor dj = 
— 2 (T;.Tj), where are generators of SU(3) corresponding to the colour charges of 
partons i and j, while the kinematic factor Wij(k) is explicitly given by the classical 
antenna function 

WW) = p ;;p s (3.H) 

where a s is defined in the bremsstrahlung scheme [33], and its argument is the invariant 
quantity k 2 ^ = 2(j>i ■ k)(pj ■ k)/(pi • pj), which is just the transverse momentum with 
respect to the dipole axis, in the dipole rest frame. Assembling these results one can write 
an expression for the average shift in the jet transverse momentum when a soft gluon 
is emitted in a selected partonic channel. These shifts must then be recombined with 
the proper weights when the full cross section is built by summing over the various hard 
scattering processes contributing to a given distribution. In any selected channel we write 

(ij) 

where we integrated over the soft gluon phase space and we defined 

5 Pt (k) = 5p+{k) e in + s P t{k) e out , (3.i3) 

with Gin being unity if the gluon is inside the jet and zero otherwise, and conversely 
for ou t- To extract the part of the soft contribution that one can associate with non- 
perturbative hadronisation effects, (5p t )h, one simply considers the region < fij in 
Eq. (3.12), with /ij an infrared factorisation scale, to be chosen so that above [i\ one may 
safely use perturbation theory. The coupling a s , which is perturbatively divergent in this 
region, is now to be replaced by a universal, non-perturbatively defined, finite quantity, 
whose moments at low energy are expected to be observable. In addition, one has to 
remove the contribution to Eq. (3.12) that would be included in fixed-order perturbative 
contributions, so that one is left with a pure hadronisation piece, which can subsequently 
be combined with fixed order perturbation theory without double counting. In Eq. (3.12) 
we replace then the coupling a s (Kt) with Sa s (Kt) = o-si^t) — o^{/k)i where a^ T (Kt) is the 
standard perturbative coupling, which can be expanded in powers of a s (pt) to the desired 
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order of perturbation theory. We note that the standard choice for [ii in the case of LEP 
event-shape studies was 2 GeV, but it should be emphasised that the sensitivity to the 
choice of this scale is only O (a™ +1 ) , if one correctly combines the non-perturbative result 
with perturbative corrections evaluated to O (a™) [13]. 

Let us continue by writing, purely for the sake of calculational convenience, 

Spt(k) = 5pi(k) + (Sp+(k) - 5pt(k)) 9 in , (3.14) 

where we used the fact that 6 ou t + @in = 1- In this way, we have divided the integral 
in Eq. (3.12) into a 'global' term, involving an integral over all of phase space of the 
unrecombined gluon contribution bp^ , and a term involving an integral over the interior 
of the jet region. We demonstrate in the Appendix that the global term does not produce 
a leading power correction, so that the complete leading contribution arises from the term 
involving the integral over the jet region. 

We note here that at the level of the single-gluon calculation all the jet algorithms 
function in an identical manner, so that a soft gluon is recombined with a hard parton 
(and 0; n = 1) if Srj 2 + 5<p 2 < R 2 . In the present case, since we have fixed the trigger jet 
at rj = (f) = 0, the gluon is in the jet for rj 2 + <fi 2 < R 2 . As noted above, the only relevant 
contribution to the present accuracy arises from the second term in the sum on the RHS 
of Eq. (3.14). As a consequence, we only need to evaluate integrals of the form 

<« t = £c«/ K d „ g *- ^ -w) _ m _ Sp;) 9iii , (w _ ^ , 

(3.15) 

for each dipole, recombined with the colour factors appropriate to the parton channel being 
considered. 



3.1 Dipoles involving the trigger jet and incoming partons 

Consider first the dipole formed by one of the incoming partons (say the one with four 
momentum p\) and the outgoing hard parton corresponding to the trigger jet, which we 
denote with the label (lj). The transverse momentum with respect to this dipole is given 
by 

2 = 2 (Pi • k)( Pj • k) = 2fe 2 (cogh ^ _ cog ^ (3 ig) 



Using 5pf(k) — Sp t (k) = k t cos<p, and changing variables from kt to Kt,ij (since the fac- 
torisation scale m sets a limit on the range of not of kt), we obtain 



<Sp,r = C„ U W ) / - , , (3.17) 

^ J v 2 +< j>2 <R 2 Ztt V2(coshr? - cos 6)2 



where we defined the first moment of the non-perturbative coupling 5a s below the factori- 
sation scale as 

1 [Vi 

A(m) = - / dK t 5a s {Kt) ■ (3.18) 

7T Jo 
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In order to make an explicit connection with event shape studies [13], we can rewrite 
Eq. (3.18) in more detail as 

A(m) = ~m 

7T 

where uq = Jq 1 a s (k t )dk t is the average coupling over the infrared region, familiar 

from event shape studies, and we have carried out the subtraction of the perturbative 
coupling, a s (pt), to two-loop accuracy in the MS scheme, where K = Ca (jg — ^Pj — 
Note that similar expressions in Ref. [13] are rescaled by the Milan factor, accounting for 
gluon decays; this rescaling will be discussed below, in Section 3.5. 

The integral over the soft gluon direction can be evaluated by choosing polar coor- 
dinates in the rj — 4> plane and expanding in powers of the radial variable. Discarding 
the spurious collinear divergence arising when the gluon is emitted along the outgoing leg, 
which cancels against an identical one in the 'global' term, as shown in the Appendix, one 
finds 

. C, AM (-1 + ±R - JLH' - ^ + O (*)) ■ (3.20) 

By symmetry, an identical result is obtained for the 2j dipole, formed by the trigger hard 
parton and the other incoming parton with momentum p2 ■ 



"o (A*/) - oi s {pt) 



ft 
2vr 



In 



El 



«a(Pt) 



(3.19) 



3.2 Dipole involving the trigger and recoil jets 

In this case the transverse momentum of the soft gluon with respect to the dipole, which 
we label with (jr), is given by 

K \ . r = k 2 t e" T ' (cosh 2 r] - cos 2 (f) , (3.21) 

which leads to the integral 

<«!f> = C jr A M [ *,» 2 ™* (3.22) 

JtP+<P<b? 27r (cosh 2 T) - cos 2 <p) 1 

This can be evaluated as before, using polar coordinates and expanding in powers of R. 
One finds 

(5p t )P = CjrAim) (-1 - \R+ -^R* - ^R 5 + O (i? 7 )) . (3.23) 

3.3 Incoming dipole 

The dipole involving the two incoming partons, labelled as (12), is the simplest to compute, 
since in this case Kt,i2 = kt- For this dipole the integration over the interior of the jet does 
not produce al/R correction, since there is no collinear enhancement for radiation emitted 
by the incoming partons as the jet becomes narrow. We expect radiation from this dipole 
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to behave essentially like an underlying event contribution, and this is in fact what we find. 
The relevant integral is 

(<*Pt)h 2) = £12*40/) / drj ^ coscp 

J v 2 +( f,2 <R 2 ZTT 

= C 12 A{^)R.h(R) 

= CuAQjli) {^-R 2 - ±R 4 + ^R 6 + O (i? 8 )) , (3.24) 
where J\ is the Bessel function of the first kind. 



3.4 Dipoles involving incoming partons and the recoil jet 

As these dipoles involve radiation which is not associated with the trigger jet, we expect, 
and verify, that they do not generate any 1/R enhancement. Considering for example the 
(lr) dipole, the relevant transverse momentum is 

K 2 t lr = 2 k\ e~ v (cosh r] + cos </») , (3.25) 

which leads to the integral 

<rf> = Or \A(„) j 4, £ .. Sm * , . (3.26) 



2 J v 2 +(t) 2 <R 2 2vr ^ ( cosh v + cos ^ 2 

This gives 

(Spt)i r = C lr AM {^R 2 + ^i? 4 - J^R 6 + O (i? 8 )) . (3.27) 

Clearly, the (2r) dipole gives an identical contribution. 
3.5 Leading power correction 

Having computed the results for each individual dipole, we can now perform the sum in 
Eq. (3.10) to obtain our estimate for the shift in the jet pt due to hadronisation effects. We 
concentrate on the leading behaviour at small R, including terms up to O(R), since at large 
values of R we expect the hadronisation component to be dominated by the underlying 
event contribution, which behaves like R 2 . Clearly, each parton channel will have a different 
Pt shift, because of the different colour weights in Eq. (3.10). Here we choose as an example, 
to illustrate our point, the scattering of non-identical quarks, 

?(pi) + q'fa) - ?(ps) + q'(p4) ■ (3.28) 

The calculation of the colour weight of the various dipoles in Eq. (3.10) is a textbook 
exercise [34]. The result is 



j^qq ^qq g 



M q Q q ^ qq 



2C F (W u + W 23 ) + — [W 12 ; W u ] ) , (3.29) 
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where, in a notation similar to [34], we have defined 



[W 12 ; W 34 ] = 2W 12 + 2W U - W13 - W u ~ W 23 - W 24 , 



(3.30) 



which is a collinear finite combination of dipoles contributing at subleading order in the 
number of colours N c . Since we are interested in just the 1/R and R terms, we need 
only consider dipoles that involve the hard parton responsible for the trigger jet, which we 
identify with the quark carrying momentum p%. Adding together the contributions from 
the relevant dipoles one easily finds 



We note several features of Eq. (3.31), most of which are common to the other parton 
channels as well. The most striking aspect is clearly the singular dependence on the jet 
radius R, arising from all dipoles involving the trigger jet. This singular behaviour is a 
reflection of the fact that, as one makes the jet narrower, one increases the pt loss due 
to hadronisation, which also explains the negative sign of the leading term. The leading 
behaviour as R — > arises from the collinear singularity of the matrix element, which is 
also responsible for the InR enhancement of the perturbative result described in Sect. 2. 
For this reason, one expects the 1/R behaviour to be accompanied by the colour charge 
of the parton initiating the jet, which in the above case is Cf, and this expectation is 
indeed confirmed when one combines the dipoles. When the trigger jet originates from a 
gluon, as it would be for instance in the gg — > gg channel, the leading 1/R behaviour of 
the hadronisation correction can thus be simply obtained with the replacement Cf — ► Ca- 
Another point to note about the 1/R term in Eq. (3.31) is the fact that its coefficient, 
2Cf-A(hi), is one half of the coefficient of the 1/Q correction to the thrust variable in 
e + e~ annihilation, as computed in [8]. In order to match more recent evaluations [35-39], 
this result for the thrust must actually be rescaled by a factor -M, where M is the Milan 
factor. Specifically, the factor ^ accounts for the use of a fully dispersive coupling (one 
expressed in terms of the gluon virtuality rather than its transverse momentum), while 
the Milan factor M accounts for the non-inclusiveness of final-state observables, i.e. the 
fact that a virtual gluon inevitably splits into massless objects, but the observable may 
change its value unless all of the gluon decay products are emitted into the appropriate 
region of phase space. For observables that are linear in multiple soft-gluon momenta, 
the Milan factor is independent of all other details of the observable, and numerically one 
finds M ~ 1.49 (for nf = 3 in the non-perturbative region), or equivalently ~ 0.95, 
which is a modest overall correction. We note, however, that jet transverse momenta are 
in general not linear in soft-gluon momenta, because of the non-trivial way in which soft- 
particles recombine with each other and into jets. The one known exception in this respect, 
among infrared and collinear safe jet algorithms, is the anti-fet algorithm 2 , introduced and 
discussed at length in [40]. Thus, while at present we can tentatively associate the size of 
our 1/R effect with the thrust calculation, this implicitly assumes that an analysis similar 

2 More generally, all algorithms with p < in Eq. (4.5). 




(3.31) 
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to that performed in Refs. [35,36,38] will yield results that do not differ significantly from 
the value computed for M in the case of linear observables. 

We finally observe that the coefficient of the term linear in R in Eq. (3.31) is numerically 
less than 20% of the coefficient of the 1/R term. This means that the 1/R behaviour should 
be a good approximation to the full result over a wide range of values of R, up to R ~ 1, 
a fact which will be useful to construct simple approximate formulae in Section 5. 

3.6 Mass distribution 

In order to illustrate the generality of the method, let us briefly summarise the results 
one obtains when performing the same computation with a different observable. We now 
compute the hadronisation correction to the jet mass, 5Mj, using the framework described 
above. Working again near threshold, the variation in the jet mass due to the emission of 
a single soft gluon is given by 

SMf (kt, V, 4>) = 2pj -k = k t Vs (cosh r) - cos (f>) , (3.32) 

provided the jet algorithm recombines the gluon with the jet. Clearly, there is no correc- 
tion if the gluon is not recombined. As the singular l/R behaviour of the hadronisation 
correction to the transverse momentum of the jet was due to unrecombined gluons becom- 
ing increasingly collinear to the originating hard parton when R — > 0, there should be 
no corresponding singular piece for the jet mass. We find in fact that the leading power 
corrections vanish at least as R for all dipoles. Specifically 

SM? 12 

5M hr 

5Mf Ar 

We observe that the small radius behaviour of the leading power correction to the jet mass 
distribution is softened with respect to the pt distribution by two powers of R. This is a 
consequence of the fact that gluons emitted inside small cones make small contributions 
to the jet mass, just as for the pt distribution. Gluons which are not recombined with the 
jet, on the other hand, in this case make no contribution, so that the singular behaviour 
at small R is absent. Note however that the contribution of the underlying event to the 
jet mass is 0(R 4 ), so that again the hadronisation component dominates by three powers 
of R, for R < 1. 

4. Comparison with Monte Carlo results 

Two simple results emerge from the analytic approach of Section 3, concerning our ex- 
pectations for the non-perturbative modification of a jet transverse momentum. First, 



- c « ^ f {* + i> 3 + ^ + jL r7 + ° (*)) • 

= CjrAUn) ^ {r+^R 5 + O (i? 9 )) , (3.33) 
= C lr AUii)^- ( -R A + — R 6 + -^-R s + O (R 10 )) . 

2 -r 25g -r 589 g 24 t- V J J 
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non-perturbative radiation associated with emission from the jet's originating parton, in a 
selected production channel, should modify the jet pt as 



{6pt)h = -Ci^A(jii) + 0(R) , 



(4.1) 



where C{ is the colour charge appropriate to the parton originating the jet 3 . Recall that, 
in the single soft gluon approximation the scale A(m) is related to the scale appearing 
in analytical studies of hadronisation in e + e~ and DIS collisions. A first test of our ap- 
proach would then be to fit data for jet distributions using expressions like Eq. (4.1), with 
a numerical value for A(/j,i) borrowed from event shape studies, including the appropriate 
Milan factor, as discussed in Sect. 3.5. We must then expect some degree of loss of uni- 
versality across different algorithms and observables, which introduce different nonlinear 
multi-gluon effects. For our illustrative purposes, however, it is sufficient to be aware that 
2Cf .A(2GeV) ~ 0.5 GeV, corresponding to the amount of non-perturbative transverse mo- 
mentum radiated per unit rapidity with respect to a qq dipole, while a gg dipole must be 
reweighted by a factor Ca/Cf = 9/4. 

Our second result is that corrections associated with the dipole of incoming partons, 
and in fact more generally with any dipole not involving the trigger jet, give a contribution 
that scales as R 2 , i.e. in proportion to the jet area. These corrections effectively mimic 
underlying event corrections, and are indistinguishable from them. Borrowing our result in 
Sect. 3.3, and denoting the scale of transverse momentum emission per unit rapidity along 
the beam direction with Ajje, we expect the jet transverse momentum to be modified as 4 



The scale Ajje, as the notation suggests, cannot be easily related to A(m), since it receives 
contributions from the underlying event, i.e. the interactions between the proton remnants. 
The functional dependence on R through the combination RJi(R) corresponds to the 
definition of 'passive vector area' of a jet given in Ref. [41,42]. 

It is of interest to compare the two sets of contributions, Eq. (4.1) and Eq. (4.2), to 
what is observed in parton-shower Monte Carlo event generators. To do so, we begin by 
generating dijet events; after parton showering, we identify jets with a chosen jet algorithm, 
and label the transverse momenta of the two hardest jets as pf^ ps , (i = 1, 2); we select events 

for which the hardest jet satisfies 55 GeV < p^ s < 70 GeV; next, we allow either the 
hadronization stage to take place, or both hadronization and underlying event generation; 
finally, we rerun the jet finder. After hadronization the two jet transverse momenta will 

(i) 

be changed, and we denote them with p t ^; similarly, after generating also the underlying 

3 As noted in the Introduction, this form of R dependence for the leading power correction was already 
suggested in [7] and its full form could probably have been deduced from [6]. It seems, however, from 
discussions of hadronisation in the literature [17, 18], that the community was not aware of the implications 
of these results. 

4 The terms of higher order in R are specific to our chosen recombination scheme, the _E-scheme or 
four-vector recombination. 



(8pt)vE = A UE RJi(R) 



Aue 
2 



(RZ-R*/8 + 0{R 6 )) . 



(4.2) 
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event, the momenta will become pJ(j E - At this point, we can define the hadronization 
contribution to the jet transverse momentum to be 

while the underlying event contribution is 

(Spt)vE = g (pIue - pIps + pJ'ue " P?ps) " {Spt)h • (4.4) 

Note that with most modern underlying event (UE) models it is not possible to carry out 
the UE generation independently from hadronization: one must first determine (6pt)h with 
the global UE switch turned off, and subsequently carry out a separate run with the switch 
turned on, in order to deduce (<5pt)uE- 

Our results are shown in Fig. 1 for four different jet algorithms, as a function of 
the jet radius R, for qq — > qq interactions (summing over the flavors and antiflavours 
of the incoming and outgoing quarks) in pp collisions, in Tevatron Run II kinematics, 
y/s = 1.96 TeV. The upper curves in each plot are for (5pt)\jE, while the lower ones show 
(5pt)h- In each case we show results from Pythia 6.4 [20], (tune A, using the default, 'old' 
shower), and from Herwig 6.5 [21] with the Jimmy UE model [22] (with the Atlas tune [43]). 
Additionally, for the hadronization results, we show the prediction based on the one-gluon 
emission approximation and universality assumption, Eq. (4.1). 

The four jet algorithms we use are all infrared and collinear (IRC) safe, as they must 
be for any analysis beyond leading order. Specifically, three of the algorithms we use 
are sequential recombination algorithms; they operate by defining interparticle and beam- 
particle distance measures via 

$ ss min(** *g) M±M , d g s fcg , (4.5) 

(v) 

and then recombining particles i and j if one of the af- is smallest, or defining particle i 

(p) 

as a jet, and removing it from the list of particles, if af^ is smallest; the procedure is then 
iterated until no particles are left. The choice p = 1 corresponds to the well-known (in- 
clusive) kt algorithm [28,29]. Choosing p = gives the Cambridge/ Aachen algorithm [30], 
probably the simplest possible IRC-safe hadron-collider jet algorithm, corresponding to the 
successive recombination of particles that are closest on the y — 4> cylinder, until all are 
separated by at least R. Its simplicity makes it our preferred choice in those figures where 
for brevity we consider only a single algorithm. Our third choice, p = — 1, dubbed 'anti-A^' 
algorithm, is novel [40], and has several interesting properties: for example, it behaves 
like a perfect cone algorithm, in that hard jets are nearly always circular, with radius R 
on the y — 4> cylinder; furthermore, it behaves linearly in soft gluon momenta, so that its 
hadronization corrections should have better universality properties, as was the case for 
event shapes in e + e~ annihilation. As a fourth example, we use a seedless stable-cone 
algorithm, the SISCone algorithm [31], with a Tevatron run II type of split-merge proce- 
dure [32]. It is similar to the midpoint algorithm [32] used at Tevatron, differing mainly in 
its seedless nature and the resulting fact that it is IRC safe at all orders. 
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Figure 1: Modification of the pt of jets due to the underlying event (upper curves) and hadroniza- 
tion (lower curves), for qq — > qq scattering at the Tevatron Run II (jap, ^/s = 1.96 TeV), comparing 
Pythia 6.412 [20] (tune A) and Herwig 6.510 [21] with Jimmy 4.3 [22]. In the case of hadronization 
the Monte Carlo outputs are compared to the analytical result, Eq. (4.1). Dijet events are selected 
containing an underlying qq — > qq scattering, and with the requirement that at parton shower level 
the hardest jet has 55 GeV < pt, ps < 70 GeV. The non-perturbative corrections shown correspond 
to the average for the two hardest jets. 



A striking feature of Fig. 1 is that all four jet algorithms have similar R dependences: 
despite much debate about the relative merits of different algorithms, they all behave 
identically in the presence of a single soft gluon, and, as a consequence, there are strong 
similarities in their non-perturbative behaviour. 

The non-perturbative contribution that is generally acknowledged as being best mod- 
elled is hadronization: indeed the Monte Carlo models agree well between each other, as 
is natural given the extensive tuning of the hadronization-related parameters at LEP. The 
analytical prediction, Eq. (4.1), reproduces both the shape and normalisation relatively 
well, though the quality of agreement varies somewhat with the algorithm: the kt (SIS- 
Cone) algorithm, for example, has slightly less (more) negative corrections than predicted. 
These differences may be a reflection of the breakdown of the single gluon approximation, 
caused by non-linearity of the kt, Cambridge/Aachen and SISCone algorithms. A full 
treatment would involve the inclusion in Eq. (4.1) of the non-trivial corrections associated 
with double-soft gluon emission, i.e. the (non-universal) Milan factor. Clearly, this would 
be of interest for further study. 

The dependence of the hadronization contribution on the jet colour factor is visible 
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Figure 2: Similar to Fig. 1, but for the gg — > gg underlying scattering channel, instead of qq — > qq. 
For brevity only the Cambridge/ Aachen result is shown. 
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Figure 3: Similar to Fig. 2, but for the LHC (pp, y/s = 14TcV) rather than the Tevatron. 

in Fig. 2, which displays results for the gg — > gg channel at Tevatron. We see that the 
hadronization contribution roughly doubles compared to qq — > qq scattering, again in 
reasonable agreement with the analytical prediction. Fig. 3, on the other hand, shows the 
independence of the hadronization contribution on collider energy, displaying results for 
gg — > gg scattering at LHC, which are easily seen to be almost identical to the Tevatron 
results. 

Underlying event contributions are much less well understood than hadronization. Not 
only is no analytical prediction available for their normalisation, but Pythia tune A and 
the Jimmy Atlas tune, despite both being tuned to similar Tevatron data, give underlying 
event correction to jet pt's that differ by a factor of two at the Tevatron, as is visible 
in Figs. 1 and 2. As can also be seen from these figures, the UE contribution is largely 
independent of the hard scattering channel (qq — * qq versus gg — > gg). It does however 
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Figure 4: Plot of the UE contribution to the jet pt, for Tevatron and LHC gg — > gg events (selected 
as in Fig. 1), rescaled by the factor RJi(R), corresponding to the passive area [41], as calculated 
for a 1-particle jet. 

depend strongly on the collider energy, as seen in Fig. 3, which illustrates the huge size 
of the UE contribution predicted by Monte Carlo models at LHC. Interestingly, though 
Jimmy and Pythia differ for the Tevatron, they agree for LHC 5 . 

Let us now turn to the functional form of the UE corrections, Eq. (4.2), and the inter- 
pretation of its normalisation. We consider (5pt)\jE divided by its predicted R dependence, 
given by RJ\(R); this should give an ^-independent result equal to Atje- The results 
are shown in Fig. 4. Whereas for Jimmy there is near perfect agreement with the scaling 
prediction, with Pythia there is clear deviation from it towards small 6 values of R, which 
can be interpreted as a correlation between underlying event activity and the properties 
of the hard scattering. Further investigation reveals that this is caused by strong colour 
reconnections in the Pythia underlying event model and tunes, an effect which is enhanced 
for gg — ► gg scattering as compared to qq — > qq scattering 7 . Note in any case that the 
physical impact of this correlation is limited, since the relative deviation from the expected 
RJi(R) form is significant only in the region where the absolute size of the UE contribution 
is modest. 

Finally, let us comment on the physical scale associated with the underlying event, Aue- 
At the Tevatron the value that emerges from the models is Atje(1-96 TeV) ~ 2 — 4GeV, 
while at LHC it is AuE(14TeV) — lOGeV. This is an order of magnitude larger than the 

5 We have also examined other Pythia tunes (DW, DWT, SO, S0A). All are nearly identical at Tevatron 
energies, while, at LHC, tunes DW and S0A are similar to tune A, and tunes SO and DWT give results 
that are about 40 — 50% higher. 

6 Note however that we are unable to consider the smallest R values, because the subtraction procedure 
in Eq. (4.4) makes it difficult to determine accurately the small UE contribution relative to the large 
hadronization effect. 

7 This feature is present to the same degree in tunes DW, DWT based on the old shower, and in all cases 
in which the old shower is used, it seems to be independent of the precise procedure used to extracts the UE 
contribution. In tunes SO and S0A, based on the new shower, it seems that the effect may be less strong, 
however it is difficult to make a firm statement because at small and moderate R the extraction of (<5pt)uE, 
eq. (4.4), is affected by non-trivial interplays between the parton shower and the UE multiple interactions. 
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corresponding scale for hadronization of a quark (2CfA — 0.5 GeV) or a gluon (2CaA — 

1.0 GeV). In all cases the scale relates to the amount of radiation per unit rapidity with 
respect to the emitter. The much larger scale associated with UE is not unexpected: one 
can either interpret it in terms of multiple gluon-gluon interactions between the proton 
remnants, or alternatively in terms of an effective perturbative saturation scale, induced 
by small-x dynamics; both would lead to substantially higher values of the effective scale 
than expected for normal hadronization. 

Though not surprising, the large value of Aue does have implications for choices of 
R in jet finding. Arguments based on purely perturbative considerations (and on nor- 
mal hadronization) suggest that R = 1 is an optimal value. The non-perturbative scale 
associated with the UE is so large, however, that for moderate pt jets it can easily be 
comparable to the perturbative contributions, proportional to a s Pt- This, for example, is 
the explanation for the poor behaviour of the kf algorithm in [44] with the 'recommended' 
R = 1 choice, while smaller R values, used more recently for the kt algorithm in [45], lead 
to more reasonable UE corrections. In contrast, when examining very high p t jets, since 
the UE should be pt independent, it will become advantageous to return to R ~ 1, so as 
to minimise the size of perturbative corrections, which behave roughly as ~ a s pt InR. 

5. Experimental considerations 

There are a number of respects in which the results presented above are relevant exper- 
imentally. On one hand, knowledge of the relative size of various non-perturbative and 
perturbative contributions as a function of R can provide guidance in the choice of optimal 
values of R. On the other hand, experimental data can provide interesting cross-checks of 
our understanding of the R dependence, for example through direct measurements of how 
the transverse momentum of a jet depends on the parameters of jet algorithm, and through 
studies of quantities such as inclusive jet cross sections. In what follows, we will sketch 
three possible avenues of experimental investigation on these issues; they should be seen as 
examples of what could be done, which will have to be adapted to the chosen observable 
and to the needs of the given experimental setup. 

5.1 Radius optimisation 

We have seen that there are three main sources of corrections to the transverse momentum 
of a jet: perturbative radiation; hadronization, i.e. non-perturbative effects associated with 
the jet itself; and the underlying event, i.e. low-p t effects associated with proton-remnant 
interactions. Each of these effects has a different dependence on the jet pt, on the colour 
factor associated with the hard parton originating the jet, on the choice of R in the jet 
finder, and finally on the collider energy, as illustrated in table 1. 

Let us examine the implications of this understanding on the choice of R in various 
experimental contexts. There are two principal scenarios to be considered, corresponding 
to two rather different usages of jets. One may use jets for the identification of underlying 
kinematic structures, such as mass peaks for the top quark or other heavy particles, or in a 
search for hadronic decays of a hypothetical Z'. Alternatively, one may wish to compare jet 
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Dependence of jet {5pt) on 






'partonic' pt 


colour factor 


R 




perturbative radiation 


~ a s (p t )p t 




hxR + 0(l) 




hadronization 






-1/R + 0(R) 




underlying event 






R 2 + (R 4 ) 





Table 1: Summary of the main physical effects that contribute to the relation between the trans- 
verse momentum of a jet and that of a parton, together with their dependence on the properties 
of the parton, the jet radius R and collider centre of mass energy. Cases labelled "-" do not 
have any dependence on the corresponding variable in a leading approximation, but may develop 
anomalous-dimension type dependences at higher orders. 



data (such as the inclusive jet spectrum) with high-order perturbative QCD calculations, 
and attempt to deduce information about fundamentals of QCD or the electroweak the- 
ory. In the first case, one seeks to extract the cleanest possible kinematic structures, and 
therefore one should minimise both perturbative and non-perturbative modifications of a 
jet pt] in the second case, one presumes that the perturbative loss is calculated with good 
precision for typical ranges of R, and one wishes to minimise the two non-perturbative 
contributions, since they cannot be precisely computed from first principles. 
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Figure 5: Contributions to the average bp\ from perturbative radiation, hadronization and under- 
lying event, for quark jets at the Tevatron, as a function of R. For the perturbative contribution 
we have used Eq. (2.7); for the hadronization contribution we have taken just the 1/R term; for 
UE we have used the full R dependence and set Aue = 4 GeV. 



When considering how well one can reconstruct kinematic structures such as mass 
peaks, one needs to know the dispersion due to both perturbative and non-perturbative 
effects, as well as any non-trivial correlations between them. Although this goes beyond 
the scope of what has been calculated in this paper, some basic quantitative information 
can nonetheless be obtained, by arguing that the dispersion on a jet pt (and therefore on 
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any kinematic structure) can be approximated by the uncorrelated sum 

(Sp 2 t ) ~ (5pt) 2 pcrt + (Spt)l + (Sptfvv . (5.1) 

This size of the summed squared contributions is shown as a function of R in Fig. 5, using 
our analytical results for the perturbative and for the two non-perturbative contributions, 
for 50 GeV quark jets, with the normalisation for the underlying event contribution set to 
Aue = 4 GeV, corresponding roughly to Pythia's estimate for the Tevatron. One observes 
that while perturbation theory prefers 8 R ~ 1, the significant underlying event contribution 
leads one to favour somewhat smaller values, R ~ 0.7 — 0.8. For this pt value, on the other 
hand, hadronization has a relatively limited effect. 

One can make plots similar to Fig. 5 for a range of jet transverse momenta, for quark 
jets and gluons jets, and at Tevatron or at LHC. In each case one can determine an 
optimal R, minimising the combination of perturbative and non-perturbative contributions, 
Eq. (5.1). The results are shown, as a function of pt, in Fig. 6. 




50 100 500 1000 

P, [GeV] 

Figure 6: The R value that minimises the sum of squared average perturbative, hadronization and 
UE contributions, as a function of pt- The approximations are the same as those in Fig. 5, except 
that for LHC we have used Aue = 10 GeV. 

One sees that in the high-pj limit, where perturbative radiation dominates, R should 
be chosen ~ 1, as expected, since this minimises | mi?|. At lower transverse momentum, 
Pt ~ 100 GeV, the contribution from the UE, with its large intrinsic scale ~ 4 — 10 GeV 
becomes of similar relative importance as the perturbative a s pt term, so that it becomes 
advantageous to decrease R, in order to take advantage of the R 2 reduction of the UE 
contribution. As a consequence, even though on average jets are slightly narrower at high 
p t than at low p t (essentially because of the reduction of a s ), it is actually advantageous to 

8 The dispersion of the perturbative contribution will not of course precisely vanish for R = 1 — rather, 
Fig. 5 should be thought of as the extra dispersion over and above whatever minimal perturbative dispersion 
would be present for R = 1. 
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quark jets gluon jets 



Tevatron 
LHC 



0.63 0.83 
0.46 0.61 



Table 2: R values that minimise the two non-perturbative contributions in various circumstances. 

use larger R values at high p t , because in absolute terms perturbative radiation is roughly 
proportional to pt- 

In contrast to this slightly counterintuitive result, when going from quark to gluon jets 
one's intuition about the correlation between jet width and choice of R is reliable: for gluon 
jets, perturbative radiation and hadronization are both larger, whereas the underlying event 
is unchanged; the optimal R, therefore, is larger. One may finally note that at LHC one 
needs a smaller R than at Tevatron, because of the noisier underlying event. 

The situation is different if one considers observables where perturbative effects are 
accounted for by higher-order calculations. One should then minimise just the sum of the 
squares of the hadronization and UE components. Ignoring all but the R 2 term of the UE 
contribution, this implies that the optimal R is given by 



where, as usual, C\ is the colour factor for the hard parton responsible for the jet. Using 
the same values for Atje as in Fig. 6, this gives the results in Table 2, which are independent 
of pt to a first approximation, and actually correspond to the zeio-pt limit of Fig. 6. 

It should be kept in mind that these results involve many approximations: for example, 
we have ignored differences between jet algorithms, we have used just the small R limit 
for some terms, we have assumed that average shifts are indicative of dispersions, and 
we have neglected the issue of actually being able to resolve separate jets for complex 
massive particle decays. We believe, nevertheless, that the basic conclusions are valid all 
the same: in reconstructing kinematic structures one should prefer larger R at high p% and 
for gluon jets, and smaller R at the LHC than at the Tevatron; when perturbative effects 
are calculated separately, on the other hand, the minimisation of the non-perturbative 
contributions alone leads one to favour somewhat lower values of R. 

5.2 Measurement of (Spt) 

Throughout this article we have discussed the change in a jet transverse momentum due to 
perturbative and non-perturbative effects. It should be emphasised, however, that this (5pt) 
is not an operationally well-defined quantity away from the threshold limit. The reason is 
that we have imagined using a parton to provide a reference transverse momentum, with 
respect to which we then discuss the shift (Spt). Beyond leading order, collinear divergences 
and quantum mechanical interference between emissions make the choice of a reference 
parton p t inherently ambiguous. Were we discussing e + e~ collisions, this would not be 
a serious issue, since the e + e~ centre of mass energy itself would provide a meaningful 
reference scale. At hadron colliders, however, and notably in dijet production, such an 
unambiguous reference scale seldom exists. 




(5.2) 
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To study directly the shift {Spt), either at NLO or experimentally, one must therefore 
find a physical way of identifying a reference scale for the jet. One possible approach, 
which we propose to follow here, is to introduce a reference jet definition, D re f, and then 
measure how a jet transverse momentum changes if one uses an alternative definition D a \ t . 
One complication is that different jet definitions lead, in principle, to different sets of jets, 
so that one must be able to identify which jet, with definition D a n, corresponds to a given 
jet with definition -D re f. An algorithm for doing so could be outlined as follows. 

1. Identify the set of reference jets that are of interest, j^J, . . . for example the 
two hardest jets in the event, using the definition -D re f. 

(i) (k) 

2. Associate with each jet in the reference set, j^, the alternate jet j^ lt with which it 
shares the most pt- In other words, find the k that maximises 

vf l) = £ Pt,e, (5.3) 

where the sum runs over particles i that are contained in both jets. 

3. Reject events in which the same alternate jet is associated with multiple reference 
jets 9 , as well as events in which a reference jet shares no pt with any of the j^, 
pJ w) = 0,Vfc. 

4. Discard alternate jets that have not been associated with any reference jets. 

Being interested, say, in the R dependence of (dpt), one might then perform the following 
analysis. Begin by selecting, with the reference algorithm, events in which the two hardest 
jets are both central, so that they are well measured (set, for example, |y^| < 2); require 
the sum of the jet pt's to be in a some limited range, for example 55 GeV < ^{p^ e {+P^ e {) < 
70GeV; repeat the jet finding with the same jet algorithm but with an alternate radius 
Rait', find the jets that match the two reference jets, and determine the difference in pt 
between the alternate and reference jet definitions, using 

For illustrative purposes, this procedure has been carried out on LHC events simulated 
with Pythia and with Herwig, and the resulting (Spt) is shown Fig. 7, both at parton level 
and at hadron level (including hadronization and UE). We used the Cambridge/ Aachen 
algorithm with reference jet radius i? re f = 0.7. One notes that the moderate differences 
between the Pythia tune A and Jimmy underlying events (cf. Fig. 3) are clearly visible 
here. One also observes that perturbative and non-perturbative effects are of comparable 
sizes. 

9 In the studies described below, we find that this occurs only very rarely: in a fraction of a percent of 
events when using the Cambridge/ Aachen algorithm both as reference and alternate, and in a couple of 
percent of events if the alternate algorithm is SISCone. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



n alt 

Figure 7: The average shift, (5p t ), as defined in Eq. (5.4)), between an alternate (Cam- 
bridge/Aachen with variable i? a it) and a reference jet definition (the same algorithm, with fixed 
-Rrcf = 0.7), determined from Monte Carlo simulation. Hadron level results include the UE contri- 
bution, whereas at parton level it has been switched off. See text for further details. 



Given a measurement of (Spt), an interesting study would be to compute the average 
p t change at NLO, {Spt)^LO, f° r example with NL0jet++ [46]. Such a calculation would 
include configurations with up to three partons in a jet (a level of detailed jet structure that 
is reached in the inclusive jet spectrum only at NNLO). Given this accurate perturbative 
knowledge of 5pt, one could then fit the two non-perturbative components, in order to get 
data-based constraints on their magnitude independent of parton shower Monte Carlos. 
Specifically, one would write 10 

(SptKRat) = (SpthtoiRsxt) - 2 (Ci) (-!- - -M A{jj.i) 

\ -Kalt -Kref / 
+ (i? a lt Jl(i4lt) - Rrci Jl(Rrei)) A UE , (5.5) 

where (Cj) is the average colour factor that takes into account the fraction of quark and 
gluon jets. Using Eq. (5.5), one would be able to fit the two non-perturbative parameters 
A(fii) and Aue- One would expect A(fii) to agree qualitatively with value extracted from 
event shape studies in e + e~ annihilation and DIS [13], and one would obtain a first direct 
estimate of Atje- 

A similar analysis could also be carried out in other contexts, for example W^+jet 
events. In that case one would consider only a single reference jet, and one could use 
MCFM [47, 48] for the perturbative calculation. Since the fraction of quark and gluon jets 
would be rather different in this case, compared to inclusive jets (especially at LHC), this 
would provide a powerful check of the dependence of hadronization on partonic colour 

10 The hadronization terms of higher order in R have been neglected here, since in the region of R where 
they contribute significantly they are dominated by the UE contribution. 
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factors, and it would also test whether the UE is independent of the underlying hard 
partonic reaction, as expected in many simple models. 



5.3 Inclusive jet spectrum 

A context where non-perturbative corrections to jets have often been discussed [17, 18] is 
the inclusive jet pt spectrum. In Ref. [17], for example, one sees that introducing a non- 
perturbative shift in pt is necessary in order to fit data for the ratio of inclusive distributions 
measured at UAl and at Tevatron. In order to estimate the effects of the nonperturbative 
shift we have computed here, let us begin by assuming that we have at our disposal a jet 
definition such that quark and gluon jets can be distinguished at parton level (an example 
of such a definition is given in [49]), and such that it can be consistently employed at 
hadron level as well (which is not the case for the definition in [49]). Then we would be 
able to write 

%, <»> - ~ <«>») + ^(» " <«>-) ■ < 5 ' 6 > 

where da t,PT jdpt is the perturbative distribution for jets of flavor i, which is then evaluated 
with an i-dependent shift 

(5pl) NP = -2 jjA(m) + R Ji (R) Aue ■ (5.7) 

It should be emphasised that, even at this stage, Eq. (5.6) relies on several approximations: 
it assumes, for example, that one can neglect the effect of event-by-event fluctuations in 
the hadronization and UE contributions, and it assumes that non-perturbative effects are 
dominated by an overall shift in the distribution. This is essentially the same approximation 
that led to the broadly successful use of a constant shift in place of a shape function when 
studying non-perturbative effects in event-shape distributions (as in [7,50,51], see however 
also [52-55]). One should note that here, even at the level of single gluon exchange, soft 
gluon resummation near threshold suggests that the distribution will be distorted, and not 
simply shifted, due to colour mixing effects [15]. It can be shown, however, that the leading 
1/R term is not affected by colour mixing in the context of NLL resummation. We can 
then safely assume that the shift in Eq. (5.7) is the dominant effect: the first term is large 
for R < 1, while the second term is large because of the large value of Aue- 

In practice, Eq. (5.6) is inconvenient because of our inability to separately define a q ' VT 
and a 9 ' PT beyond leading order (LO), for jet algorithms that are valid also at hadron 
level. We may however further simplify Eq. (5.6) by expanding in the small relative shifts 
{5pl) NP /pt- To this end, let us define 

dlnjda^/dpt) _ da^ I da PT 

dlnpt dp t J dp t 

which represent respectively the power-law fall-off of the distribution, which behaves as 
p^ n ' , and the fraction of jets that are of species i. We can the write 

dj t m = -*r {pt) I 1 + n «~ir + °[<-ir)) + ( q ~ 9 )> (5 - 9) 
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We note that the neglected term in Eq. (5.9) is expected to be of the same magnitude as 
terms that would arise from the replacement of the shift with a shape function. Finally we 
rewrite Eq. (5.9) as 

dpt 1 dp t V * f q n q + f g n g \ ' Pt ) ) ' 

In Eq. (5.10), one can evaluate da PT /dpt up to NLO (or beyond) since it is summed 
over jet flavors, while for the purpose of calculating the (small) non-perturbative shift, 
fi and rii can be evaluated at LO, where the jet flavour is unambiguous. Generalising 
Eq. (5.10) to include not only 1/R hadronization terms, but also subleading powers of R, is 
a nontrivial task: one would have to sum not over q and g jets, but rather over all different 
scattering channels (qq — ► qq, qg — > qg, gg — > 55, etc.), and furthermore the threshold 
approximation suggests that colour correlations would make it necessary to go beyond the 
shift approximation. These subleading terms would therefore be quite interesting from 
a theoretical viewpoint, while they are unlikely to have a significant phenomenological 
impact, given the likely experimental constraints in the difficult environment of hadron- 
hadron collisions. For phenomenological applications, on the other hand, a further and final 
approximation can be made on Eq. (5.10), if one assumes that quark and gluon species 
have a similar power-law fall-off, n q ~ n g . One can then simply write 

ip~ t {pt) ^i^( Pt ~ {5pt)Np ) ' {6pr) = fq {6p9t)Np + h {6p9t)Np ' (5 - n) 

Again, this is can be extended, if one wishes, to separately consider all 2 — > 2 scattering 
channels, however a full treatment of the O (R) terms in the hadronization correction will 
probably require that the prediction be rephrased in the language of shape functions. 

To conclude, we observe that in addition to the inclusive jet spectrum directly, one 
may also examine the relative difference between the distributions that emerge when using 
two different jet definitions, in analogy to what was done in Section 5.2. One may define 

P (D, a ,D M , Pt) ). (5.12) 

Perturbatively, with current tools [46], the numerator can be studied up to a,, because the 
o? s term vanishes; in other words, this is effectively a three-jet observable, which can be cal- 
culated to NLO, and thus this is a quantity that can be studied to higher perturbative order 
than the inclusive jet spectrum itself (for which a\ would correspond to NNLO). Eq. (5.12) 
is similar in this respect to the difference in p t between alternate and reference definitions 
discussed in Section 5.2. The separation of p into perturbative and non-perturbative parts 
reads 

p{D rcf ,D alt ,p t ) = p PT {D Tcf ,D alt ,p t ) (5.13) 

+ — f q n q ((5pt) N p,alt - (<frt>NP,ref) + fg n g ((<*P? ) N p,alt ~ (<fr? >NP,ref) 

Clearly, these techniques can also be applied to ratios of distributions measured with dif- 
ferent centre-of-mass energies (as in [17]), or differing with respect to other parameters. In 
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this way one may attempt to minimise or isolate the perturbative contribution, and get a 
better handle on our non-perturbative inputs. When these are established, one may apply 
the resulting shift to the distribution itself, in order to optimise the analysis in view of 
other QCD or new physics studies. 

6. Conclusions 

We have performed a first broad analysis of non-perturbative QCD effects on jet observables 
at hadron colliders such as Tevatron or the LHC, examining hadronization effects and 
underlying event corrections, and studying how they compare to perturbative contributions. 
We have used both analytic models, within the framework of the dispersive approach, and 
different Monte Carlo tools, implementing different IRC-safe jet algorithms, belonging to 
both the 'cone' and 'sequential recombination' families of jet finders. Finally, we have 
considered different observables, focusing mostly on the shift in the transverse momentum 
of a jet due to QCD radiation, but giving results also for the jet mass and for the jet pt 
distribution. 

Our most significant results are displayed in Eq. (4.1) and Eq. (4.2), and collected in 
Table 1. They show that the transverse momentum of a jet of radius R is modified by 
hadronization effects in a manner proportional to 1/R, while underlying event corrections 
grow like R 2 , increasing with the jet area. A similar pattern holds for our analysis of the 
jet mass, where in general the correction is suppressed by two powers of R with respect to 
the pt analysis, but with the same hierarchy between hadronization and underlying event. 
As we briefly reviewed, the perturbative i?-dependence for the jet pt is logarithmic at small 
R. We thus have three different functional forms for the .R-dependence of the three main 
contributions to the jet momentum. This leads to the possibility of selectively minimising 
chosen combinations of these three effects, by a suitable choice of R, in order to optimise 
a given physics analysis. 

Our results are formulated in terms of two different non-perturbative parameters: the 
first, A(p,i), describes hadronization, and we expect it to be closely related to similar 
parameters measured in event-shape studies in e + e~ annihilation, a conjecture which is 
consistent with the Monte Carlo results; the second parameter, Atje, describes the flow 
of transverse momentum per unit rapidity due to underlying hadronic activity in a given 
collider environment: its value must be extracted from data and is expected to be uni- 
versal across different observables in a given collider, while it should scale with a power 
of the centre-of-mass energy. A further characterisation is provided by the dependence of 
transverse momentum on the colour charge of the parton originating the jet: in our ap- 
proximations, both perturbative and hadronization contributions are proportional to the 
colour charge, while the underlying event is independent of it. 11 . 

These result are rather strikingly confirmed by our studies with Monte Carlo simula- 
tions. Hadronization corrections, as given by non-perturbative models embedded in Herwig 
and Pythia, indeed grow at small R, broadly agreeing with the analytic estimate both in 

11 Work is already in progress [42] to understand how perturbative and UE effects interact, so as to 
introduce a subleading dependence on the jet colour factor in the UE contribution. 
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shape and normalisation, while the corresponding underlying event corrections (computed 
with Jimmy for the Herwig simulation) grow with R as expected. 

Finally, we have given three examples of experimental studies where our results could 
be either tested or used in order to perform physics analyses. Specifically, we have outlined 
how our results are likely to impact the single-inclusive jet pt distribution, a signature 
observable at hadron colliders in general and at LHC in particular; we have indicated how 
the pt shift from parton to hadron level could be experimentally accessed in order to test 
our calculation; and we have suggested methods and choices in order to optimise the value 
of the jet radius in view of different experimental situations. 

We regard these suggestions as just basic examples of the studies that could be per- 
formed. More generally, we believe that our study emphasises the necessity for variety 
and flexibility in the experimental choices of jet definitions and parameters. In the early 
days of QCD, when theoretical predictions were mostly at leading order in perturbation 
theory, it was reasonable to adopt a rough definition of a jet, and several such definitions 
proved sufficient to gather striking evidence for many basic features of QCD and of the 
standard model. Those days are past: QCD is now precision physics, with NLO calcu- 
lations forming the standard benchmark for predictions, selected observables having been 
computed at higher orders, and all-order resummations and power correction studies also 
available. From an experimental point of view, the challenges of the high-luminosity, high- 
energy environment of the LHC will require that all the experience acquired with previous 
accelerators, and all our theoretical knowledge, be put to good use. 

We believe that in order to do that it will be necessary to take advantage of the 
full flexibility offered by modern jet tools. There are now several well-tested IRC-safe jet 
algorithms, both of cone type and based on sequential recombination. They depend on 
parameters, notably the radius R, which provide the experimenter with handles to control 
the flow and impact of perturbative and non-perturbative QCD radiation. Exploiting the 
adaptability of these jet tools will be crucial for model testing, for the validation of exper- 
imental procedures, and, last but not least, for studies that will further our understanding 
of the strong interactions. 
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Appendix 



Here we compute the global term involving the integral of the unrecombined gluon contri- 
bution 6p^~ , defined in Eq. (3.9), over all of phase space, to show that after accounting for 
the evolution of parton distributions it does not contribute any linear power correction. To 
illustrate the point, we consider the dipole formed by the two incoming legs, discussed in 
Section 3.3. To restrict our attention to soft emission alone we impose a cutoff, requiring 
that the longitudinal fraction (with respect to the beam) of the emitted gluon momentum 
k never exceed e c , implying \r]\ < r/ max (fc t ) = In e c yfs/k t . The relevant integral is then 

(<^>gi 2) = ~~ dk t 5a s (k t )dr]— (cosh(r/) + cos(f) . (1) 

Here we have set the colour factor of the dipole C\2 to be 2Cf, since we are interested only 
in the radiation collinear to the incoming legs, which we take to be quarks. In actual fact 
the colour charge Cp for each leg will be built up from a sum of dipole contributions, all of 
which involve the incoming leg in question and give the same result in the collinear limit: 
the replacement of the dipole colour factor by 2C F anticipates this fact. To illustrate our 
point, but without loss of generality, we use a toy model for the non-perturbative coupling, 
setting 5a s (kt) = AS (kt — A). In this model the quantity A, which governs our leading 
power corrections computed in Section 3, is just proportional to A. After performing the 
integral over the full range of azimuth (f>, the global term gives 

(12) c F r p™* 

($Pt) 'gi = / 6a s (k t )dk t / dr] coshry 

2C F f 

= — / 5a s (k t ) dk t smhr) max (2) 

= --L. €c y/t + 0[-= 

As we shall see, the dependence on the cutoff e c will be removed by a similar term due to 
the evolution of the parton distributions. Note then the absence of any term of order A, 
and hence the absence of a leading power correction due to this global term. 

In order to verify the factorisability of the collinear divergence in Eq. (2), let us now 
turn to the contribution of parton distributions. In the calculations presented in section 3 
we worked in the threshold limit, effectively approximating parton distributions with their 
bare expression at threshold, qo(x) = 6(1 — x). In our non-perturbative approximation, the 
distribution evolved to a hard scale Q is then given by 

q(x, Q 2 ) = 6(1 - x) + ^ f 1 [ Q ^ 5a s (k t ) (6(x - z) - 6(1 - x)) , 

Jl-€c 1 ~ Z JQo fc * 

= S (l-x) + ^ f 1 J?-(6(x-z)-8(l-x)), (3) 

TT Jl-e c 1 - Z 

where Qq is an arbitrary reference scale, we have placed the same cutoff as in Eq. (2) on the 
maximum emitted longitudinal momentum, and we have used the infrared limit (z — ► 1) 
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of the quark splitting function P qq (z). Note that when one considers parton evolution 
one automatically includes collinear branchings which push the hard scattering away from 
threshold. 

Taking into account parton evolution, we are then forced to include a further non- 
perturbative shift in the transverse momentum of the jet, given by 



dxidx 2 (q(xi,Q 2 ) q(x 2 , Q 2 ) - qo(xi) qo(x 2 )) 

__ 2C F ( f 1 dzx (v^T- 



_ w > +(1-2) (4) 

TT \Jl- ec 1 - Zl 2 J 

= -^e c VS + 0(e 2 c ) . 

7T 

One observes that the collinear divergence in Eq. (2) has precisely the form required to 
be absorbed in a 'renormalised' parton distribution, such as the one computed in Eq. (4). 
Eq. (2) then assures us that the remaining finite terms do not contain leading power 
corrections. 

Similar calculations can be carried out for all dipoles, and in each case one observes 
that no 0(A) term arises from 'global' integrations over all phase space. For each dipole, 
furthermore, collinear divergences can either be factorised into the parton distribution, as 
done above, or they cancel, as in the case of outgoing legs, against an identical divergence 
in the 'in-jet' contribution, as expected from the collinear safety of our observables. 
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